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Blue CDM-photon isocurvature perturbations are attractive in terms of observability and 
may be typical from the perspective of generic mass relations in supergravity. We present 
and apply three theorems useful for blue isocurvature perturbations arising from linear spec¬ 
tator scalar fields. In the process, we give a more precise formula for the blue spectrum asso¬ 
ciated with the axion model of 0904.3800, which can in a parametric corner give a factor of 
0(10) correction. We explain how a conserved current associated with Peccei-Quinn sym¬ 
metry plays a crucial role and explicitly plot several example spectra including the breaks 
in the spectra. We also resolve a little puzzle arising from a naive multiplication of isocur¬ 
vature expression that sheds light on the gravitational imprint of the adiabatic perturbations 
on the fields responsible for blue isocurvature fluctuations. 
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1. INTRODUCTION 

Single sealar field inflationary models generate approximately adiabatie, seale-invariant, and 
Gaussian primordial density perturbations Il2]-[T0]|. This is eonsistent with the Cosmie Mierowave 
Baekground (CMB) measurements [[TT] - [22]| and the Large Seale Strueture (LSS) observations 
Il23l [24]| . However, non-thermal eold dark matter (CDM) seenarios sueh as axions []25] - [27]l and 
WIMPZILLAs |[28] - [33]l naturally have observable CDM-photon isoeurvature perturbations (e.g. 
Il34l - l49ll f sinee the CDM never thermalizes with the photons. Indeed, it is remarkable that a sub¬ 
dominant dark matter eomponent as small as 10^^ of the total dark matter eontent ean leave an 
experimentally deteetable effeet through eosmology (see e.g. I[50l f. Lurthermore, isoeurvature 
perturbations are interesting since it can generate rich density perturbation phenomenology. Lor 
example, unlike standard single field inflationary scenarios, degrees of freedom responsible for 
isoeurvature perturbations are able to generate large primordial local non-Gaussianities [[50] - [78l . 


3 


Scale-invariant isocurvature spectrum is well constrained as its power on CMB length scales 
has to be less than about 3% of the adiabatic power [dH US l79] - [84ll . Because the largest scale 
invariant isocurvature effects on CMB measurements occurs on long length scales (e.g. see the 
appendix in [l85l l. one expects scale invariant isocurvature effects to be well hidden in any future 
observations probing short length scales. However, if the isocurvature spectrum is very blue, then 
isocurvature effects that are hidden on long length scales may become large effects on short length 
scales (see e.g. [[Ml - (89l ). If such strongly blue spectral index isocurvature signal is uncovered in 
the future, one may ask what one will learn regarding the high energy physics of the isocurvature 
sector. 

One answer to that is given by |II| in which a supersymmetric axion model is constructed giving 
rise to a blue spectrum. In that work, the phenomenologically relevant axion isocurvature pertur¬ 
bation amplitude da is assumed to be given by the frozen value dajcp-^ at horizon crossing where 
(p^ is the classical value of the radial field that breaks the Peccei-Quinn (PQ) symmetry during 
inflation. Unlike in the conventional minimal axion scenario in which the order parameter (p^ is 
sitting at its potential minimum during inflation, its value is initially displaced from the minimum 
and is slowly decreasing towards the stable minimum during inflation. Hence, the assumed frozen 
value at horizon crossing increases for larger wave vector k modes which leave the horizon later. 
In that way, a blue spectrum is generated over a k range that depends on the spectral index which 
controls the amount of time 9 + takes to settle to its minimum. Furthermore, because supergravity 
structure generically induces a Hubble scale mass [|90ll for (p^, the spectral index can be easily 
extremely blue. For example, it has been claimed that this scenario allows an isocurvature spectral 
index of n = 4 for k G kmax ] specified by kmax/kj^jn ~ exp (10). 

Partly motivated by this result, we formulate three elementary “theorems” regarding isocurva¬ 
ture perturbations with a very blue spectrum for non-thermal dark matter fields such as the axions 
with a displaced (possibly time dependent) vacuum expectation value during inflation.^ Theorem 
1 defines a superhorizon conserved quantity for systems possessing an approximate symmetry 
associated with linearly perturbed system. The merit of this theorem compared to previous discus¬ 
sions of this topic in the literature (see e.g. [|^[92ll l is its ability to go beyond the end of inflation 
and the reheating process. Theorem 2 describes under what averaging conditions that fluid quan¬ 
tities behave as dXnadIXQ- This second theorem merely restates what is known in the literature 

* The proofs are only at the rigor of a typical physics literature. 
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in the context of the theorem. Theorem 3 describes the computation of the quantum isocurvature 
perturbations. The merit of theorem 3 compared to the previous discussion in the literature is the 
explicit canonical quantization in the presence of linearized gravitational constraints. Furthermore, 
we point out the clear conditions under which the simple analytic estimates are valid. 

We then give couple of applications of our theorems. First, we improve on the naive quanti¬ 
zation of axions in scenario of BB and compute 0(\-\- ) corrections to the spectrum of with 

spectral index n. In the process, an interesting application of conserved PQ symmetry current is 
made, which explains how two independent dynamical degrees of freedom behave as a single one 
during a finite time duration of interest. The theorems also help to set a precise boundary of where 
the simple analytic computations are invalid. For example, contrary to claims of [|TI|, ^ = 4 spec¬ 
trum cannot be generated in their scenario. Another consequence of understanding the boundary 
is that if the ratio of axion isocurvature blue power amplitude to the adiabatic power amplitude 
is at most of the order of a few percent on the largest observable scales and can be described in 
terms of quantization methods presented in this paper (and implicitly approximated in [B), most 
of the cold dark matter must be made of different species. We also illustrate through example 
plots, phenomenologically interesting parametric corners of the model (having a six dimensional 
parameter space). Although observable spectra can contain breaks, these break regions typically 
contain k-space domains for which the simple analytic computation is invalid. We identify how 
large the expansion rate H during inflation can be in this class of models generating a large blue 
spectrum. A measurement of tensor-to-scalar ratio at the level of r = 0(10“^) will disfavor this 
class of models, at least in its simplest form. 

In another application, theorem 3 is used to explain why the isocurvature blue spectrum does 
not have a simple lower bound suggested by a naive operator product analysis. More explicitly, 
the isocurvature perturbations are defined to be a contrast of the form ^ C\dx — C25(^ where 
Ci are background field dependent coefficients and 50 is the inflaton field and 5;t; is the field 
responsible for the existence of isocurvature perturbations. In other words, the isocurvature field 
is always dressed with the inflaton sector. The quantum correlator {SS) would then naively have a 
piece that is proportional to C|(5050) coming from the dressing. This piece for a blue spectrum 
is of order of the adiabatic perturbation power spectrum. If the cross correlation piece does not 
precisely cancel this piece, {S^S^) would be of the order of adiabatic spectrum, leading to a simple 
lower bound. However, the theorem shows that the power spectrum of generically behaves 
independently of the adiabatic spectrum. The more broad lesson encapsulated by theorem 3 is that 
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the gravitational coupling of dx to makes dx grow an inhomogeneity that looks like 50 such 
that becomes independent of 50. 

The order of the presentation will be as follows. In Sec. three simple theorems and couple 
of corollaries useful for spectator dark matter isocurvature spectra are presented. In Sec. a 
couple of applications of the theorems are given. One application corresponds to improving and 
elucidating the computation of BH. The second application corresponds to understanding how 
dressing effects coming from the definition of the isocurvature perturbations do not mix inflaton 
field quantum fluctuations with the dark matter field quantum fluctuations because of the secular 
growth imprinting an adiabatic inhomogeneity to the dark matter field. We close with a summary 
and thoughts on future work to be done in this direction. In the appendix, we collect some results 
useful for the theorems. 

2. USEFUL SIMPLE THEOREMS FOR BLUE ISOCURVATURE MODELS WITH A SLOWLY 

ROLLING TIME DEPENDENT VEV 

2.1. Definitions 

In this subsection, we define the language used for our theorems. 

Metric and Fourier Conventions Although theorems that we present are gauge invariant, we 
will have the occasion to use several gauges in our proofs. The Newtonian gauge scalar perturba¬ 
tions will be parameterized as 

ds^ = {\+m>^^'^)dt'^-a^{t){\F2^^^^)\dx\^. ( 1 ) 

We consider slow-roll inflaton field 9 scenarios in which superhorizon adiabatic perturbations are 
approximately conserved. Conserved adiabatic curvature perturbations on superhorizon scales is 
given in Newtonian gauge by the solution [fTOl l93] - (951 

(L k) - = k = constant (2) 

where H = d/a and \k/a\ H and we have introduced the Fourier convention 

Q{t,k) = J d^xe-^'^^-^Q{t,x). (3) 

In the Newtonian gauge, the expansion is manifestly isotropic. Furthermore, has the intuitive 
interpretation of being the gravitational potential in the Poisson equation. Because of these prop- 
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erties, the field equations in the Newtonian gauge are convenient to work with when working with 
classical equations. 

On the contrary, the spatially flat gauge is more useful for quantization during inflation (see 
e.g. [l96l k The scalar metric perturbation convention in spatially flat gauge can be chosen to be 

ds^ = (1 -\-adiF^^f^dtdx^ —a^{t)\dx\^. (4) 


As shown in Sec. 2.4 the relevant interaction action derived from solving the gravitational con¬ 
straints in this gauge consist only of local terms of the fields unlike for the corresponding equations 
in the Newtonian gauge. Thus, the quantization of fields and the investigation of the subhorizon 
mode functions are technically simpler in this gauge. Hence, we will employ the spatially flat 
gauge only for quantization during inflation which establishes the initial conditions for the late 
time classical equations. 

Linear Spectator Isocurvature Field Let linear spectator isocurvature field be defined as a 
canonically normalized scalar field % = XQ{t) -f 5x^^\t,x) for which 


5p 


f o<SzW + 0(5z'"'>-) 


^0 (^during inflation) ^ ^during inflation 


dp 


(N) 

X 


0 


^^dominant 


dT 


{N)0 


< 1 


(5) 

( 6 ) 

(7) 


dominant o 

in Newtonian gauge where the subscript “dominant” corresponds to the energy density component 
that dominates For example, during inflation, “dominant” corresponds to the label (p while 

during radiation domination, “dominant” corresponds to the label 7 representing the relativistic 
degrees of freedom. Because we will focus on 


1 


W = -mh\ 


( 8 ) 


Eq. 0 translates to 


H^Mp 


(9) 


where e is the inflationary slow-roll parameter. If the inflaton potential is given as V(p{(p), then 

2 


£ = 


IV ip 


^<p(*Po) / 


( 10 ) 


where (p = (po{t) + d(p. The effective expansion parameters are Xo/Mp and slow-roll parameters 
of the inflaton field if m/// ~ 0 ( 1 ). 
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Spectral Conventions The gauge invariant speetrum of linear speetator ;t^-photon isoeurvature 
perturbations useful for Boltzmann equations is often defined during radiation dominated universe 
through 


“ 2;r2 


r d^k' 


{2ny 


{dsJt,k)dsJt,k')) 


^sy = 3(Cz “ Cy) 


r _ ^(A^) Ibaekground smoothed ^ ^Py^^ 

'2 / /A I D \ ^ 


( 11 ) 
( 12 ) 

3(Pz+'fz)time ’ 3(py + Py) 

where P,- are pressure quantities eorresponding to —T‘- eomponents of the energy momentum ten¬ 
sor. Here, the “time” average in the denominator of the definition of eorresponds to a time 
average over m"^ time seale. The “baekground smoothed” in the numerator of the definition of 
eorresponds to averaging over time seale all quadratie terms in the baekground Zq(^) appear¬ 
ing in the numerator. The variables and are eonserved outside the horizon if the pressure of 
the eonstituent is a funetion only of its energy density. In partieular, ^y eorresponds to the gauge- 
invariant eurvature perturbation if we assume that radiation behaves as a single eomponent fluid 
eoming from the inflaton deeay. For single-field inflation, observational normalization of 


A|(ko = 0.05Mpe“^) = A^Jko) 


V<p(ko) 

24n:^M^£ 


2.4 X 10^ 


(14) 


eorresponds to the eurrently known approximate value of adiabatie eurvature perturbation ampli¬ 
tude. 

The ;t^-photon isoeurvature speetrum often eontain k-spaee domains whieh ean be parameter¬ 
ized as 




n—1 


(15) 


where n is the speetral index. The isoeurvature speetrum is blue when n > 1. The primary foeus of 
this paper is regarding speetra for whieh n — 1 > 0(0.1) whieh beeome parametrieally insensitive 
to the inflationary slow roll parameter values of 0(e) < 0.02. As far as the phenomenologieal 
bounds are eoneerned, note that 


A^(k) = co^Al(k) (16) 


where (O;^ < 1 is the fraetion of eold dark matter that is in the X field as is explained in Appendix 
B The eurrent phenomenologieal bounds on A^(k)/A^^(k) for seale invariant power speetrum is 
approximately a few pereent [fT3l[TH[79] - [84ll . 
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With these definitions and assumptions, we ean eonstruet a useful statement that ean be used to 
set elassieal equation boundary eonditions before Eq. Q breaks down. The most important of the 
three theorems that will be presented below is theorem three. Note that one of the key merits of 
the theorem that we are presenting is its applieability eonneeting eomputations during inflation to 
variables during radiation domination. 


2.2. Theorem 1: Classically Conserved Isocurvature Quantity 


Here is a statement of the first theorem. In slow-roll inflationary seenarios, the linear speetator 
isoeurvature quantity 

= (17) 

where 

^Xnad = (18) 

on superhorizon length seales is approximately eonserved as long as x interaetion is dominated 
by = n^-x^ jl and gravity, anisotropie stress effeets ean be negleeted, and attraetor behavior 
of ^Xnad and Xoi^) is relevant during inflation (i.e. non-pathologieal boundary eonditions are 
ehosen for the homogeneous field) with an expansion rate of H. A suffleient eondition for attraetor 
behavior with non-pathologieal boundary eondition is 

IZo(^iniital)l ~ ^o(^initial)/^ 


vNk > 1 


(19) 


where 


3 / 4 

2 V 


( 20 ) 


and Nk is the number of efolds between the time of k-mode horizon exit and the end of inflation. 
Here, we have defined 


dxif(l,k) = j ‘Itait) + fdoxoit) (21) 

where = 0 in the Newtonian gauge (i.e. G = N) and in any other gauge G is related to the 
Newtonian gauge eoordinates through -f ((^®, Furthermore, is a gauge 

invariant quantity. Note that we have introdueed a faetor of 2 in the definition of 5;^ for later 
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convenience. Finally, note that this theorem is formulated at the classical solution level. The error 
in the conservation coming from the assumption of attractor behavior can be estimated as 

attractor error ~ O (exp [— 2v/4]) (22) 


where the coefficient of the error depends on details of initial conditions of both the homogeneous 
mode and the perturbation mode at the beginning of inflation. The fractional error O(^) in the 
conservation is approximated to be the terms that are dropped in making this statement: 


= exp [— 2v/4] + ■ 






(23) 


^^dominant 

We also implicitly assume that the post-inflationary cosmological history consists of smoothly 
connected patches of power-laws. 

It is important to note that this theorem makes conserved independently of the details not 
stated in the theorem, including some of the details of the end of the inflation, reheating, early 
radiation domination, and how Xo makes the transition from a slow-roll field to a coherently os¬ 
cillating one. In particular, the classical conservation here is valid even when e —)■ 1 at the end 
of inflation, unlike the spatially flat gauge quantity Sx^^^^Xo which undergoes generically under¬ 
goes time evolution at the end of inflation. The conditions stated in the theorem can be understood 
as a decoupling limit of the isocurvature perturbations, and this theorem establishes a classically 
conserved quantity in that limit. Note one of the important points for this paper: the numerator 5x 
and the denominator Xo must correspond to the same dynamical degree of freedom that responds 
to the same potential dominated by the mass term. Finally, note that when we state the as¬ 
sumption that the mass term and gravity dominate the interactions, we are stating that perturbative 
interactions are too weak to thermali / e the system. 


proof 

Consider the equation of motion for the perturbation variable in the Newtonian gauge 5x^^'^ in the 
long wavelength limit in which we can neglect the gradient terms: 

§X^^'^ + 3H5x^^^ +V”{Xo)5x^^^ + 2V;^{Xo)^^^^ = 0. (24) 

Here, we have assumed that gravitational interactions and potential self-interactions V^{Xo) dom¬ 
inate the interactions. If anisotropic stress effects can be neglected, the ij component of Einstein 
equations imply 




( 25 ) 
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The 00 component of Einstein equation in Newtonian gauge partially determining is 




IMj L 


^Px ^‘’Pdominant 


where 




-xl^ + XQdt5x + V'5x 


( 26 ) 


(27) 


and ^P^ominant ^ ^^dominant o dominant contribution to the energy-momentum tensor as 


discussed in Eq. In the limit 




^ (^dominant 


< 1 , 


(28) 


we see that is independently of 5x- Hence, with the condition of Eq. (281, the dependent 
terms in Eq. ( |2^ are external sources terms. Due to dilatation diffeomorphism gauge solution that 
lifts to physical solutions [|97ll . there exists an adiabatic solution 

2 : 0(0 


SxiT = -Ci 


a{t) 


dta{t) 


(29) 


where is the usual time independent gauge-invariant curvature perturbation constant determined 


by the inflaton sector cp approximately independently of dx^^'^ as long as Eq. (28) is satisfied. 
Erom the perspective of the classical equations we are discussing here, is simply a constant 
parameterizing a solution to Eq. (24) where the gravitational potential is given by Eq. ([^.^ 

Given that Eq. ( [24] ) is a second order differential equation, the most general solution corre¬ 
sponds to two independent solutions hi 2 to the homogeneous equation added to the particular 


solution given by Eq. (29): 

^ +C2h^2^ + Sxia^ (30) 

where ci and C 2 are coefficients independent of time. Hence, we see that in the limit that k/ (aH) —)■ 
0 can be neglected, the numerator of 


^Xnad 

Xo{t) 




ci{k)h^^\t,k) + C2{k)h2'' ’ {t,k) 
Xoit) 


{N)> 


ad 


Xo{t) 


(31) 


^ These statements can easily be covariantized, but such formalizations tend to obscure the intuition rather than 
to illuminate the intuition. Since our aim is to illuminate the intuition of the simple physics, we will leave the 
presentation in the explicitly gauge dependent form. 

^ Although we have not made any explicit assumptions about the background energy density, Eq. (28 1 does depend 


on the background energy density. 
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is governed by the same equation as the denominator if = n?"x}I2\ i.e. 


Henee, we ean also write 


hf'^ + = 0 

+ m^Xo^'^ = 0 - 


Xo = eihf^^ +e2h^^\ 


(32) 

(33) 

(34) 


We know that one mode deeays faster than the other during inflation. This is what we usually 


eall the attraetor behavior during inflation [|951 . We will eall the less deeaying mode h 
quantitatively, in the dS approximation, we have 


(N) 


h 


{N) 


h 


(N) 


= e 


2Hvt 


More 


(35) 


(36) 


In this ease, we thus have at the end of inflation (when this relative growth ends) 

Xo{t) ei+e2(k)0{e-^^^k) 

whieh is independent of time Nk (the number of seale faetor efolds between k mode horizon exit 
and ) as long as 

viVfc>l. (37) 

Henee, the error in the eonservation eoming from the attraetor assumption is O (exp [—2v/4]) • 
Finally, under the gauge transformation + ((^°, we have 


Sx‘'°\uk) = Sz‘"H<7)+?“aoZo(!) 




W/ 


(38) 

(39) 


on long wavelengths. This means 




ad 


(40) 


Xo{t) Xo{t) 

for general gauges G non-singularly eonneeted to the Newtonian gauge N. We thus see that this 
quantity is gauge invariant. 

Note that one may wonder whether there are other interaetions besides mass interaetions that 
would lead to the same result. To see that this is not generieally possible with only potential 

^ad 

V'^{5x)^V';{5x)5x (41) 


modifieations, note that for — ^xid^ to behave similarly as Xo{t), a generie eondition is 
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which can easily be solved to obtain 

Vxidz)^Cidx^ + C2 (42) 

which means that x interaction is dominated by the mass term. 

A trivial corollary of this theorem is to discuss the situation when the constant m is replaced 
by m{t) which is constant during a finite time interval during inflation and makes a transition to 
another value during inflation. 

corollary 1 In the context of theorem 1, suppose m is not a constant but makes a transition to 
another value during inflation: 


m^{t) 


{ mi t <tc 

m2 t >tc 


(43) 


where the transition time region near t is assumed to be much smaller in time than ^ 
and H ^. The quantity is still conserved as long as sufficient time has passed during the 
t<tc period to be in the attractor approximation just as in theorem 1: i.e. 


li{tc-tk)v{mi) > 1 


(44) 



(45) 


(46) 


again with the neglect of any possible secular effects that depend on unusual cosmological 
histories. 


proof 

On superhorizon scales, we have just as in theorem 1 proof 

= 0 (47) 




( 48 ) 
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Attractor behavior during t <tc gives for the solutions 

dXnad ^ Cl/lf)(l+O(e-2''(™0/^(4-4))) 


(49) 

(50) 


in the language of the proof of theorem 1. Using the well known “sudden” approximation, one can 
match these solutions valid for t < to those valid for t > tc'. 


e\ 


Cl 



I > \ / El] 

f h'"* <> \ / Cl \ 

Ur)<)),_ icJ 


(51) 


(52) 


where Ef and Q are independent solution coefficients specifying the Xo ^nd SXnad solutions (re¬ 
spectively) in the time region t > tc and are independent solutions in the t > tc time region. 
Clearly, we have in the region t > tc 


^Xnad _ ^^1 

Zo 


(53) 


which is a constant. 


2.3. Theorem 2: Gauge Invariant Isocurvature Spectrum During Radiation Domination 


Here is a statement of theorem 2. The radiation dominated period linear spectator isocurvature 
perturbation spectrum defined by Eq. ( [TT| ) on superhorizon length scales is given by 

cfik' 






-{S^{t,k)S^{t,k )) 


(54) 


which is time independent as long as X interaction is dominated by = nrx^ 1'^ gravity, 

anisotropic stress effects can be neglected, slow-roll attractor behavior of the Zo(0 during infla¬ 
tion is relevant (e.g. boundary conditions close to slow-roll are chosen for the homogeneous field), 
and m 3> 3H/2 during the radiation dominated time period when one wishes to evaluate this ex¬ 
pression. An important part of the linear spectator requirement is given by Eq. Q. An interesting 


point of this theorem is that is defined in Eq. (11) with Xq in the denominator while Eq. (54) 
is proportional to Xq in the denominator. 


Proof 
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Consider the eomputation of 

during radiation domination when m S> 3H/2. Sinee Xoi^) is eoherently oseillating, its energy 
density consists mostly of non-relativistic energy density. Hence, the definition of 

SXnad = Sx^^^ - Sxif (56) 


gives 


= 


Zo^f^Zfiadlbackground smoothed 5"^^'!«^lbackground smoothed 


^^o)time 

according to Eq. ( [T3] ) where we also defined the term “background smoothed.” 

(G) 


(57) 


From the definition of the conserved quantity Si in Eq. (17), we can make the substitution 




(58) 


Hence, the isocurvature classical quantity during radiation domination simplifies to 

1 + 


8 


"^^(^o)ti 


time 

^^o)time 


(59) 


where we have used theorem 1 in keeping Sy constant, and the time average is given by 


(P +^)time 


(^o)time 
^ (^o)time- 


(60) 

(61) 


This gives 


5.,=S;t. (62) 

We thus conclude that {k) is asymptotically time invariant during this radiation dominated time 
period when the conditions of theorem 1 are valid approximations. 


2.4. Theorem 3: Quantum Correlator of Linear Spectator Isocurvature Perturbations 

In addition to the conditions of theorem 2, if Bunch-Davies boundary conditions to the infla¬ 
tionary quantization are imposed, m < 37/(during infiation)/2, and a slow-roll inflationary phase 
characterized by the slow-roll function £ ^ 1 defined hy H ~ —eH^ during inflation occurs, the 
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spectrum of linear spectator isocurvature perturbations on superhorizon scales during radiation 
domination is given by 


~4 


(63) 


Tt ) \2TtxQ{tk^)) Uoy 

where tk satisfies k/a{tk) = H(tk) during inflation, gQ{m/H) is an order unity function which 
vanishes when m/// —)■ 0, and 

v = , (64) 




with ko being a fiducial wave vector (typically taken on CMB length scales). Here, this expression 
is only valid when 


viV^>|ln/| 


(65) 


where Nk ^ {f/^k, H{te — tk)}. Here, tk is the time of the horizon exit of mode k, tg is the end 
of inflation, and / is the fractional error tolerance in the computation. Moreover, for v to be real 
and also for the correction to 2^ |r(v) | to be within the error tolerance throughout the time period 
between when the mode k > ko leaves the horizon and the mode ko leaves the horizon, we must 
have 


2£fcn 


m 






< / for k > ko- 


The bounds on v sets the largest valid isocurvature spectral index of 


( 66 ) 


R = 4 — 2v 


(67) 


coming from Eq. (63). Finally, since Zo{tk) was approximated by expanding about tk^, we must 
have 


/ 


ko exp (- <k <ko exp ( — 

^ko J V £fco 


(&■ 


( 68 ) 


The fractional error 0{S) in Eq. (63) is approximated to be the terms that are dropped in making 
this statement: 


S = exp [-2v{mi)Nk] + 


5P, 


+ + ^ + (69) 


dominant 


3v^//2 Mp 


Mn 


Because we are neglecting O^Eko) contributions, this result is applicable to very blue spectra in 
which 3 — 2v S> £yto, and because go(0) = this result is applicable to 3 — 2v = 0 case. 
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Before moving onto the proof, let us make some eomments. One immediate implication of this 
theorem is that large isocurvature spectra with blue spectral indices are easy to generate for masses 
of order H. For example, the supergravity 7]-problem for the inflaton can be translated to the 
isocurvature sector to generically expect rr?/H^{tkQ) ~ 0(1) during inflation. On the other hand, 
this does not mean that such spectra are measurable since any observable must be multiplied by 
the background energy density of %, which tends to dilute away when m/H is not small (as noted 
by [I98ll l. Still, there are examples [[U of supergravity models generating measurable, strongly blue 


spectrum. This will be discussed in depth in Sec. 3.1 In many applications, a/c 1 which means 


that Eq. (691 imposes a constraint on |;to(^fc) I for tho longest wavelengths leaving the horizon. 


proof 

To make an inflationary prediction based on theorem 2, we need to match the classical gauge 
invariant quantity 

C (fik' 4 ^ ^ 

J (70) 

to a quantum correlator computation. A technical difficulty lies in the fact that dXiiadi^ik) is 
not the field that is being matched to the quantized field (recall it is only a particular part of a 
classical solution whose initial condition statistics are captured by the correlator of Eq. (fTO])). One 
way to quantize is to quantize the modes of the linearized equation of motion for 5x^^\t,k) in 
a particular gauge G. It is important to note that we only need the correlator value computed 
from field quantization during the few efolds of the horizon exit since by that time the correlator 
evolves as a classical statistical object and (k) is frozen according to theorem 1. If this were 
not true, we need to understand the full time evolution of the quantum correlator which is a messy 
task. Below, we will quantize in the spatially flat gauge (G = sf) instead of the Newtonian gauge 
(G = N) because the gravitational potential does not appear in the coupled oscillator equations. 


After quantizing, one can compute the correlator of Eq. (70) by computing 


Eet us now consider the details obtaining Eq. (70) by matching to the quantum computation 
during a few efolds of horizon exit. Since the relationship between the spatially flat gauge isocur¬ 
vature field and the Newtonian gauge field is 




( 71 ) 
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we have 


Xo Xo 





> (I) 


, -29t 

.H J Xo \ Xo 


nyZ, 

Aq 

P 


,W^W\ 
k p I 


(72) 


Beeause of the linear spectator property, during inflation is given by the adiabatic mode. We 
parameterize ^X^f^ using the nonadiabatic mode ^X^nad^^^P)'- 

Sx^P - SxP{t,p) + dxPi^,p) = -Cp^ jdta + 8xpl{t.p) 


4> 




Cp 


H r 

1 -/ dta 

a I 


(73) 

(74) 


We find 


' ^Xp^ Sxp^ 

Xo Xo 


^xPdi^,k) 5xPdi^,p)\ I ^ \ 
Zo Zo / + 

^xPit,l) 


' Xo " 

HXo, 


-291 


^nad ' 

Xo 


^HXo, 


(75) 


To extract pxladipk)dxPit,p)'j using the spatially flat gauge correlators, we also need an 
expression for pxiad(^^^)^p} terms of pxp^■ Following a similar procedure as before. 


this is given by 


" " 'szri(<7)fe)-(Ci&> 



Hence, we have arrived at the desired expression 


Xo 

HXo 


(76) 


^ pxPdi^,k) SxP{t,p) 
Xo Xo 



(77) 


where the left hand side is interpreted as a correlator of stochastic classical fluctuations and the 
right hand side is computed at the quantized level. Note that this expression indicates what is 
conserved is effectively the naively expected quantity IXo with the gravitational potential 


effect subtracted. One may naively think that the second term of Eq. (77) contributes significantly 
when m/i7 ~ 0(1). However, we will show below how this term is exactly canceled out. 
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Let us now compute the correlators explicitly in the leading slow-roll approximation. We recall 
that in the spatially flat gauge, the mode equations are homogeneous with off-diagonal mixing 






(78) 


which allows one to use adiabatic approximation to quantize the system since the mass matrix 
varies sufficiently slowly in time during slow-roll period for approximate diagonalization. Mak¬ 
ing a time independent rotation to the approximately diagonal basis, we find the approximately 
diagonal modes hyj satisfy 






= 0 


(79) 


where 


-I- ^ (h^ — — ^0 


+ 2 


v;(Po 

MjH 


<pixo 

2M4//2 


{3r\y -6£)H^ 


(80) 

(81) 
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^<PoXo ^ (PoXo ^ (poXo 


Ml 






-2^-2 


= V" -I- ^ — — ■^0 

2M4//2 


+ 2 


and U^M^U is a diagonal matrix. Let’s define 


K = 


Ml, 


?=o 


MjH 


Kxo 

MjH 


^Mp 
<piXo 

2M4//2 


XoVlp 

Mp 


PeXo 


mm+0(§. 


M„ 


sgn((Po) 


(82) 


(83) 


(84) 


L = 


Mk: 

Mf, 




|f=0 


(85) 


(3r]v-6£)H2 

where k I because of Eq. 0 , As far as L is concerned, consider two cases. One case will be 
\L\ S> 1 and the other case will be |L| <C 1. The former is the more interesting case, since the other 


case is the one which people have phenomenologically considered most commonly and has been 
well established. 
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Case with |L| S> 1 


^2 0 2 


+ 


•M 


U22=l 




+ 


K-213 

L 


lljc^ 3 k ^ 


+ A, 


U 


12 


= -kA„+5M+x3 


•M 


U 21 = kAm- 


A 


'M 


•M 

3k^ k 

~Y~1? 

3k^ k 

~Y~1? 

f2 


8 

11^4 


8 


2 L 2 

3^2 


+ ... 


+ ... 


1 +Am 


9k‘2 


2 L 

-A^( 

' K 

9^3 ^ 

L3 

2 L , 


+ ... 


+ ... 


( 86 ) 

(87) 

( 88 ) 

(89) 


with formal perturbation power Am assignment as Mf^{0) = (9(Am), = 0{Xm), ^ 22 ^^) ^ 

0 ( 1 ) and t = 0 is defined to be initial time here when the modes are deep within the horizon (i.e. 
when k/a{0) 3 > // ). Note that this diagonalization differs in a numerically insignificant manner 
from that used in [|99ll where the diagonalization is carried out the horizon exit time (See Eq. ( 8 ) 
in [|99l ).^ Hence, because 

Vn 

(90) 


^«i. 

Mp 


the off-diagonal element of the mixing matrix U is certainly generically negligible at the initial 
time. 

Based on the smallness of mixing in U, one might naively think that the mixing continues to be 
negligible for the far superhorizon evolution of if Xq/^p 1 since 


+{^+Ml^] Sx'-’f'’ = 0(Zo/Mp)5<pWI 


(91) 


which may naively allow us to neglect the dependent right hand side in the limit Xq/^p 0 . 


However, this is not quite right. Since K ^ 1 and L S> 1 which implies 


(92) 


^ Since the vacuum physically corresponds to zero on-shell particle states when the modes are approximately 
Minkowskian, diagonalization of the mass matrix when the modes are deep within the horizon is physically more 
faithful. On the other hand, as long as the mass matrix is sufficiently slow in its time dependence, the distinction 
is not numerically important. One might also argue that a time dependent diagonalization in the spirit of WKB 
approximation is even more appropriate when the modes are deep within the horizon because as shown in OlOOI . 
the time dependent rotation always generates mixing at any time jo/<Po ^ 0. However, as noticed in Il99l , the ratio 
of the off-diagonal mass squared and the diagonal mass squared is proportional to slow-roll parameters whose time 
variation is suppressed by higher order in slow-roll. In that sense, neglecting the time dependence of the mixing 
matrix is justihed. 



















20 


the field will eventually grow relative to to see this, use the trial superhorizon 

positive frequeney mode solution e~'^* and consider the resulting dispersion relationships in 
the constant H limit. In contrast, in the decoupling limit, will evolve independently of 

5;^^Because of these two features, the mode obtains an effectively inhomogeneous 

(i.e. sourced) contribution proportional to 


'''particular^ ' 


(Po(f) 


(93) 


where we have used slow-roll equations of motion but have not used the superhorizon approxima¬ 
tion.^ The fact that this is a valid solution in the subhorizon region is explicitly demonstrated in 
Appendix]^ On top of this, one can add a homogeneous solution 




(94) 


where h 


(• 5 /) 

z 


is a free oscillator field satisfying 


2,M) 


d^h 


dt^ 


+ +( 3 + 3^12 1 = 0 


'-X 


with the normalization 


{hx^\t,k)h^x^\t,p)) d^^\k + p) 




(95) 


(96) 


and C is a coefficient determined by Bunch-Davies boundary conditions (note the basis h defined 
here already satisfies the Bunch-Davies boundary conditions). Note that if H and M 22 are con- 

(sf) 

stants, hy ’ is composed of Hankel functions. 

Another way to justify the decomposition into the particular and homogeneous solution here is 


that we have to keep the second term of Eq. (94) which goes as 






—m 


Zo 




(97) 


<Po(0 3y/2£H^sgn(poMp 

even in the decoupling limit of Xo/^p ^ 1 because of the unsuppression due to the relative grow¬ 
ing nature of relative to ■ This contribution is small on subhorizon scales but grows 


^ Neglecting 5^6/) influence on 5^6^^ is certainly valid since the mixing is small and 5^6'/) grows relative to 
5^6/) when Eq. (92 1 is satisfied as one can see by considering the positive frequency mode equations. 

® In spatially flat gauge, this particular solution is related to the gauge invariant perturbations through 1 ^ « 
/(po. 
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after the modes leave the horizon. This is a type of seeular effect in which long time behavior 
unsuppresses the naively suppressed contribution. 

Let’s now estimate C in Eq. (|9^. By definition 




(98) 


where and are decoupled at the initial time t = 0, with mass matrix eigenvalues m\ 


and respectively (see Eq. (jv^). To determine C, we want to match Eqs. (94) and (981. Note that 
the decoupling used in Eq. ( [^ already drops contributions except for those that can 

eventually be unsuppressed by the relative growth of compared to after the modes 

leave the horizon.^ In the subhorizon region, there is no enhancement due to relative growth of 
to ■ Since we will be matching in the subhorizon region, we can drop the U 12 and 

U 21 terms when matching: 




At earlier times when the modes are subhorizon, from the normalization ( [96| ) we have 


if 


kja-^ 


ml ^ MI2 


(99) 

( 100 ) 

( 101 ) 

( 102 ) 


consistently with Eq. ( [92| ). Then Eq. ( [99| ) gives 

\C\^ ~ (5Xd^\t,k)dXa^\t,p) 


+ 


2:0(1) Y 




(■?/)/ 


(103) 


Assuming {^8Xj^\t,k)dXj'^\t,p)'^ ~ (^d^j'^\t,k)d^jf'^\t,p)) in the kja —)■ 0 ° limit, we can 
conclude that 


C= 1 + 0 


Zo 

3V^H^Mp)' 


(104) 


^ Although 1121541 ^ in Eq. (98 1 is also unsuppressed by 5(Pj^k^ in the superhorizon region, the coefficient U 21 = 


0(y/eK) has an extra power of slow roll parameter e when compared to Eq. (97 1 . We do not really need this fact 


for the demonstration here since we are matching in the subhorizon region. It is being mentioned to note that the 


second term of Eq. (94i is not coming from U 21 which is something that is being evaluated only at the initial time. 














Eq. thus can be written as 


Sz 






Zo 


where we used that in spatially flat gauge 


Zo 

HZo 




Zo 


( 105 ) 




(106) 


Hence, approximating the Bunch-Davies state correlator near horizon exit time as 

2 






a-^ 




d^^\k + p) 


where 


^ , 4 

^ ~ 2\l 


(107) 


(108) 


we expand this in the limit k/ (aH) —0 in the usual manner to arrive at 


Zo Zo 

which is only valid when 


‘ 2’-l-»|r(v)|^(A)"'^| (109) 




»|ln/| 


( 110 ) 


where / is the fractional accuracy desired and H is approximately constant. Even Eq. (110) is 
about justifying keeping the leading term in the Hankel function expansion, we will soon discuss 
that this also corresponds to the attractor behavior needed for the validity of theorem 1 that will be 
used here. 

Eq. (1091 is not manifestly frozen since for example the first term seems to dilute as a^. How¬ 


ever, just as we stated in theorem 1 (which is more general for classical solutions than this com¬ 
putation of a few efolds during inflation), this expression does freeze. To see this, during the few 
efolds period of horizon exit during inflation, we can solve 


Zo + 3//Zo + m Zo = 0 


( 111 ) 


to obtain the attractor solution 


a{t) 


— ^+V+O{£)go{m/H) 


Zo(i) ~ Zo(?fc) 


( 112 ) 
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where k/a{t]^) = H{tk), is a funetion whieh vanishes in the limit m/H —0, and the approx¬ 
imation drops the deeaying solution. Henee, one sees that all the a(t) dependeneies eaneel in 


Eq. (1091 and are replaeed by a{tk) ■ 

Let us now eonsider matehing the quantum eomputation to the elassieal variable of theorem 
1. With the usual Buneh-Davies normalizations, 5Xnad on superhorizon seales is determined by 
purely imaginary mode funetions: i.e. the time dependenee of mode funetions are determined by 
by Hy^\k/{aH)) where 

(113) 

where the —)■ 0) —)■ oo. Thus 5Xnad eommute like a elassieal random variables onee 7^^ 

dominates over 7^^ The time seale on whieh 7^^ dominates over 4^^ is 


X = 


IvH 


(114) 


This is eoineidentally the same time seale as the elassieal attraetor behavior diseussed in theorem 


1 and approximately the same time seale after whieh Eq. (110) beeomes satisfied. Onee 


exp[-{t-tk)/x] </ 


(115) 


eondition is satisfied, one ean mateh the eorrelator of elassieal random variable to the quantum 
eorrelator on superhorizon seales using Eq. ( |77| ): 

(5A) = 4(2.)3sW(? + p){^2^(>'-»|r(v)|^(±)"''^}. (116) 

In other words, to use theorems 1 and 2 to eonneet late time isoeurvature to the quantum eorrelator. 


we must satisfy Eq. (110). 


Einally, it is instruetive to eonvert Eq. (116) into the form of Eq. (54). We start by using 
Eq. ( |112| ) to reeast Eq. ( |116| ) as 


A4(fc)=4 


(117) 


. xi{tk){2n)^ 

where one must keep in mind that despite its appearanee, we have already evaluated this at a time 


when a{t) S> a{tk) and we have already assumed Eq. (110) is satisfied (we will find the eonstraint 


imposed by the assumption below). Choose a time eorresponding to a fidueial mode horizon 
exit time to write 

H{tk)=H{tk,){^\ ( 118 ) 

















24 


implying 


A4«=4 


a{tk) 

22v//2(^,J|r(v)|2 


-e-l 


^ \ (3—2v)—2e+0(e)go(m//f) 
ko 


(119) 


( 120 ) 


Xi{tko){27t)^ 

Note also that we are keeping the ln(A:/ko) enhaneed e dependenee while dropping other non- 
enhaneed e dependenees. Sinee we only speeify that gQ{m/H) vanishes when m = 0, the —2e 
power is numerieally meaningful in the eurrent estimate only when m = 0, whieh is really not 
about the blue speetra. Nonetheless, we keep it here to eonneet this speetra to the massless axion 
speetra. 

Let us now find the parametrie region where one expeets the attraetor to be reaehed eonsistently 


with Eq. (110). An obstaele to satisfying Eq. (1101 in the eontext of Eq. (1171 is that either inflation 
ends too early or the assumption that H is eonstant during the realization of Eq. ( |115 1 is invalid. 
The assumption of eonstant H breaks down on a time seale At satisfying 


\H\At = fH 


( 121 ) 


where / is the aeeuraey desired. During slow-roll, we have 


At = 


eH' 


( 122 ) 


Henee, Eq. (110) beeomes 


vA,»|ln/| 


(123) 


where A4 G {//C/t, H{te — h)} where ti^ is the time of the horizon exit of mode k and tg is the end 
of inflation. 

Note also even if \H\At/H is within error toleranee /, it may still be bigger than what is need 
to keep V real. Expressed in terms of v with H evaluated at the fidueial ko horizon erossing, we 
need 




(124) 


whieh translates to 

trP' 

2 £fco ^2('^^ ^y 2 ^^ <1 for k > ko. (125) 

Eurthermore, the prefaetor terms 2^ |r(v) | and Xoih) were approximated by expanding about 
tkQ assuming a eonstant H. We require the eorreetion to 2'^ |r(v) p to be within the error toleranee. 
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This means that in considering the large enhancement situation 2'' |r(v) p ~ 1/v^ near v = 0, we 
impose 


2''Mr(v,)|2^4( l + 2e, 


m 


-ko 




Inik/ko) ) < :^( 1 +/), 


(126) 


that gives Eq. ( 661 . Finally, Eq. ( |122 ) restricts k through the horizon crossing considerations to be 
in the range 


/ 


kQ exp (- <k <ko exp — 

^ko J V 


(a- 


(127) 


Case with |L| ^ 1 In this case, we have 


,=_y|2 

"n 




f=0 


( 3 t ] v /- 6 e )//2 


< 1 . 


(128) 


Because we still have Xo/^p ^ 1 ^nd 


M; 
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rn^XoV^ 

Mn 


sgn^o, 


we have the hierarchy 




(129) 


(130) 


Hence, the off-diagonal element of the mixing matrix U is still suppressed. (In fact, it is now sup¬ 
pressed relative to all of the matrix elements.) Furthermore, again because of the mass hierarchy in 


Eq. (130), does not grow relative to unlike the situation explained just after Eq. (911. 


Hence, we can in this case completely ignore the mixing. 


Finally, we see the term shifting the field in Eq. (77) is negligible to leading slow-roll 

order accuracy: i.e. 


Xo 




■m^\ 

CV I - I 


37/2 ; 


< (3T]K-6e)Cfe. 


(131) 


Hence, we again conclude that the isocurvature spectrum is given by Eq. (117). 

Just as for theorem 1, a trivial corollary of this theorem is to discuss the situation when the 
constant m is replaced by m(t) which is constant during a finite time interval during inflation and 
makes a transition to another value during inflation. 

corollary 2 In the context of theorem 3 and just as in corollary 1, suppose m is not a constant but 
makes a transition to another value during inflation: 


n?{t) = 


mi t <tc 


m2 t >tc 


(132) 
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where the transition time region near t = tc is assumed to be much smaller in time than 
and H Suppose sufficient time has passed during the t <tc period to be in the attractor 
approximation just as in corollary 1: i.e. 

iV(t„t^)v(mi)>|ln/| (133) 


where N{tc,tk) e {//Cyt, {tc-tk)H}, 


4 m? 




(134) 


and / < 1 is the error tolerance in the computation. The spectrum is still given by Eq. (63) 
with V —)■ v(mi) for modes k in the range 


^min < ^ < ^max 


(135) 


where ^max is the smallest k that among [k that saturates the inequality of Eq. (133), 
^oexp(//e^J,fcoexp } and 


^mm = ) ^oexp ( 


f 


(136) 


where t/, is the beginning of inflation. The fractional error 0{S') in Eq. (63) receives contri¬ 
butions from 

\Xoitk)\ , \Xoitk)\ 


^ = exp[-2v{mi){tc-tk)H] + 


SpP 


+ 


5pW . ^ ' 3y^//2 Mp 

'^dominant 


+ 




(137) 


The error contribution proportional to 1 /a/c in Eq. (137) can be rewritten as abound on Xoitkp,)/M^ 


if we require it to be less than / and require kn^ ~ uqUoTz/I^ mode to be the quantizable within 
the current analytic treatment: 




M„ 1 






|Zo(4o)l 


(138) 


This and Eq. ( [T^ give 






)2v-l 


|r(v)p 




Tt 


2 nr 

3flP 


(139) 


kHg is the wave number corresponding to the observable universe today that in a comoving box with a length of 
L = 4/flo^^o- 
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Since we know that the left hand side is eonstrained to a few pereent level (at least for the seale 
invariant speetra), phenomenology requires that the dark matter fraetion in x be very small for a 
very blue speetra: i.e. (0;^ <C 1 for m/// n-~j 0(1). 

proof 


For those modes whieh ean have Buneh-Davies boundary eonditions at early time and exits the 


horizon long before t = tc will lead to a speetrum where Eqs. (109) and (110) are valid with 
v(m) —>■ v(mi): 

« (lnfiV.)Ck + p) [ inv(^i))l^ 


Zo 




n% 


aH 


I 


VW- 2 YI 9^2- 


If the dominanee of over 7v(mi)(-^) oeeurs as 

exp[-(t-tfe)2v(mi)i7] -C/, 


(140) 

(141) 

(142) 


with t<tc, elassieal behavior is heuristieally justified, and one ean mateh the quantum eomputation 


to the elassieal solution. Sinee the eonstant H approximation result of Eq. (140) requires t < 


tk + f /(SkH) and Eq. (142) must oeeur with t < tc, we arrive at the eonelusion of Eq. (133). Also, 
for V to be real throughout the time period between when the mode k> ko leaves the horizon and 
the mode ko leaves the horizon, we must have 

m 


v(mi) > sJ^£k^{Nk-Nk^) 


//(ko) 


for k > ko. 


(143) 


whieh follows from the justifieation of Eq. (125) with v —>■ v(mi). This sets an upper bound of k 
to be at 


k < ko exp 


2£^a 




(144) 


The k bound eoming from x{h) being eonneeted to xiho) through a de Sitter spaee solution is the 
same as Eq. 


3. APPLICATIONS 

3.1. Improvement of the Axion Blue Isocurvature Scenario |IT]| 

In this seetion we apply our theorems to the seenario of [HI and eompute 0{\/{n — 4)^) eorree- 
tions to their blue speetrum. 
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3.1.1. A Review of the Axion Blue Isocurvature Scenario mi 


We begin by reviewing [[T]|. They eonsider a renormalizable superpotential 

W = h(4>+4>_ 

where the subseripts on indicate t/(l) global charges. The F-term potential is 

Vf =I (I P + I P) I ^0 P • 

A flat directions of Vf exists along 

4>+d> = d>o = 0. 

Their soft-SUSY breaking terms are assumed to be 

Vscft =m^|4>+|^ + m2 |d>_|^ + mo|4>oP 

where m,- = O(TeV). The Kaehler potential induced potential is 

y/f = C+//2|$+|2+c_//2|^_P + CqH^ I d>0 P 


(145) 


(146) 


(147) 


(148) 


(149) 


where c+,_p are positive 0(1) constants. In addition to these, there can be FI induced trilinear 
terms which can spoil the flat direction. Hence, they assume that the inflaton sector can be arranged 
to have H Fa such that the flat directions are only lifted by the quadratic terms. 

Looking along the flat direction of Eq. ( lAl) (more explicitly, setting d>o = 0), they have the 
effective potential being 


y ^ I 12 + C+//21 14>_ 12. 


(150) 


During inflation, the minimum of d>± lies at 




— ) Fa. 


(151) 


They assume d>± starts out away from the minimum with a magnitude larger than this and ap¬ 
proaches the minimum during inflation. This implies the f/(l) symmetry is broken during infla¬ 
tion. Hence, there will be a linear combination of the phases of which will be the Nambu- 
Goldstone boson associated with the broken [/(I). Hence, they make a judicious sigma model 
parameterization 

(?>± [ . 

d>± = — exp i - 


72 "V72«P± 


( 152 ) 
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where (p± and a± are real. For our explanation later, keep in mind that (p± and a± are four distinet 
dynamieal degrees of freedom. 

The potential in the new variable is 


V ^ cos 


a+cp- +a-(p+ 

V2(p+(p- 


+ h^F^F-h^cp^_cpl 


+ -c+H^(pl_ + -C-H^(p^ 


(153) 


For any fixed value of (p±, one ean see that a linear eombination of a± will be the Nambu- 
Goldstone boson with an approximately flat direetion. That axion eombination ealled a is de¬ 
termined in their equation 15: 


a = 


<P+ 


(pl + cp^ 


-a+ — 






(154) 


whieh is not to be eonfused with the seale faetor (whenever it is not elear from the eontext, we will 
add the subseript fiigcale 1° denote the metrie seale faetor). The orthogonal eombination 


b = 


9 


(pi+ 9^ 


za-\- -|- 


9+ 


has a potential 


Vt, = —h^Flcp^cp^ cos 


9l + 9^ 


' (pi+ (p^ 


(155) 


9+9- 


(156) 


2FI along the flat 


(157) 


They assume is initially large (of order of Mp). Note that sinee (p+(p- 
direetion, if (p+{ti) ~ Mp, then 

Henee, the field b during this time near ti has a mass of order h(pp hMp 3> H (for Mp) and 
thus is assumed to be settled to the minimum of Z? = 0 during inflation. The mass squared matrix 
of (p± also says that9+ mass squared is of order of sueh that it ean be dynamieal. Sinee b is 


sitting at Z? = 0, Eq. (155) implies 


5a ~ ——Sa-L. 
9+ 


(158) 


Henee, the angles appearing in Eq. (152) are 


56. = 


5a4 


V 29 + 


(159) 
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and 


de = 


da 


V2(p- 


Because of Eq. (157), we have the fluctuations in the a field being 


a ~ a+ 


= 0+a/2(P+. 


(160) 


(161) 


They therefore claim that the quantum fluctuations of 


(162) 
a 

is frozen upon horizon departure. They assume a is massless and assign 5a an amplitude of 
H / (2;r) where H is approximately constant during inflation. With that reasoning, they write 


H 


^amplitude 


(163) 


It is important to note here that as far as identifying da with H/ {In) is concerned, they are ne¬ 
glecting the kinetic term induced mass of the a field which is what we are going to correct below. 
Note (jo+decreases as a function of time as 


"“"’•“'(Si) 


+v 


with 


2 V 9 7/2 


(164) 


(165) 


where ^ is the mass squared of the (p^ field and 0 being a fiducial time. Since horizon exit 
condition gives 

^scale^^i) _ 


which gives 


^amplitude 


^scale^^) \^scale^^i)^^^i) 


H{h) 


-e-l 


(|-v)(l+e)- 


-(^1)0+ V^scale(^i)^(^i) - 
This is indeed a blue spectrum in the limit v 0 (i.e. for m/H —)■ 3/2). 


(166) 


(167) 


3.1.2. Improvement 


According to our theorem 1, what needs to be frozen is da/a. Since a is decaying with time 
(in contrast to the standard axion scenario), da cannot be massless. On the other hand, [ill assigns 
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a massless spectral amplitude to 5a of H/{2n) (see Eq. (163) above). The goal of this subsection 
is to remedy this. 

Although the action is a bit simpler if we choose other c-model parameterization, we will 
consider the explicit axion parameterization {a,b, (p±} considered in [[l|| for transparency in con¬ 
necting with this work. The kinetic term term is 


K = |(94)+|2 + |,9d>_ 




|2 

\2 


(168) 


-dfia [sin^ yd^ In (p_ -|- cos^ yd^ In (jo+j -|- bdp,a [d^ In cp —d^ In cp^] sin ycos y 


— ^djj^b [cos^ yd^ In(p_ -f sin^ yd^ In(p+j -f 


d^(p-d^(p- 


1 


I 2-2 ab rL . 

-H-^cos rsin y -^cos rsinr-f 

[2 A(pI ^ ^ (pI ^ ^ A(p 

1 


/ 2 . • 2 \ COS'^ 7 

—j (cos 73-4 sm 7 )- 


sin ^7 


+ 


-dfx(p-d^(p+^ [—2a(?cos'^7-fa^cos^7sin7—3(?^cos^7sin7-l-2a(?cos^7sin^7] -|- 


2 ' ■ ■ <P+ 

d^(p+d^(p 


1 4 S ^ • 2 .. 2 \ • 2 

-H-^ COS 73 — ^cos 7sm73-^ (sm 7-|-4cos 7 sm 7 

2 A(py (py 4(pj^ ^ 


(169) 


where 

cos7= - ■ (170) 

^ 9++ (P^ 

Because b and cp have masses much larger than the expansion rate H, they will sit at the minimum. 
Although the global minimum of this potential is at 


b = 0, (p+\,ni„ = V2^- = (171) 

assuming Fa is sufficiently larger than H, when cp-^ is displaced from its global minimum, the cp 
will sit at the cp variation minimum of 


(P 


2E3 


1 


<P+ 13-2 


£_ 172 

(p2 


(172) 


which varies as a function of time because (p^ varies as a function of time. On the other hand, 
because its mass is heavy, its fluctuations are not dynamically important. 

For later discussion of the dynamics of the (p+ and a, it is useful to note that the a field here 
corresponds to coset space parameterization of the spontaneously broken PQ symmetry: U{1)pq 
defined as 


d>± 


(173) 
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Because U{\)pq breaking VEV is considered in the dynamics (i.e. (p± dynamics), the full dynam¬ 
ics is characterized by an unbroken U(l)pQ. U(l)pQ gives rise to an exactly conserved current 
JpQ that is conveniently expressed in terms of a±: 

JpQ = + -a^d^cp-. (174) 


This current is one piece of crucial information not discussed in [|T1 . 


Integrating out b (i.e. setb = 0 in Eq. (169)), we have 


1 


K -(^a)^ —[sin^7^^1n(j!)_-l-cos^7(9^1n(p+]-I- 


da(p-d^(P- 


1 Cp' 2 . 2 

- H-^ COS 7sm 7 

2 4(p2 ‘ \ 


+ 


^df^(p-d^(p+^ [a^cos^7sin7] -h 
1 


2 

d^(p+df^(p 


- H-^ cos 7 

2 4(p2 \ 


(175) 


Next, we will focus on the dynamics while 9 + S> <P- and 9 + 3> a where we will be able to integrate 
out (p . Since 7 1 in this regime, we can then drop the sin 7 terms to write 


K 


1 






1 a 

- + 


2 n 


.2 4(p2 


da(p+d^(p+. 


(176) 


Next, since the mass of cp is of the order 


h^(pl > H^, 


(177) 


we can integrate out (p^, leaving a and (p+ as only dynamical degrees of freedom. Note also that 


from the consideration of the kinetic term Eq. (176), the axion is no longer shift invariant as long 


as p 0. Note also that the simplification embodied in Eq. (176) breaks down when 9 + is 
near its global minimum of Eq. ( 171| ). Also, in order to stabilize (p+ at the minimum, we have 
implicitly assumed 

, /c2IcZ ^ 

(178) 


f:- > 




and in order to decouple 9 and b at the minimum, we have assumed 



\/c~c+ 2 

h'^{c-+c+) 


(179) 
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The equation of motion for a and are 


D(p+ + ^d, 


Vs ^ I2(p 




^ □o)+ 

□a — a -= 0 

(P+ 

1 1 


[y/gad^a] +d(p^V = 0 


(180) 


(181) 


2 <?>+ ^/g 

Next, consider the background and linearized equations for a and (p+ with the cp and b sitting at 
its minimum 

b = 0 


(P- = 


2F^ 

<p+ 


1 




(182) 

(183) 


leading to the approximate potential of 


VI 


f 2 Ac-h^F^ 

“ T V+’’+ A2^|+2c_H2 


(184) 


With this truncation, we have the conserved U{1)pq charge density being 




PQ 


4^4 

'— f- + 1 ) ( 9 +d+ — a+(p+). 

n 


(185) 


Because of the spacetime expansion, this charge density dilutes away: i.e. 

^+(0 _ <P+(l) 

«+( 0 ) ^ <?>+(o) 

to exponential accuracy. Note that during the early period of inflation when (j!)+ ^ Fa, we have 


(186) 


(187) 


a ~ a+. Eq. (186) and theorem 1 explains why during this time period 

da <P+(ti) da 
—- ———— = constant. 

a{t) a+{ti) q)+{t) 

In other words, da/ (p+ is frozen because it behaves like da/a because of U{1)pq even though the 
numerator and the denominator are independent degrees of freedom for which theorem 1 would 


not always apply. For the complete spectrum calculation, it is better to rewrite Eq. ( |186[ ) in terms 
of a{t): 

+■?+(') 


a{t) 


(188) 


Eet’s consider the magnitude of the initial charge that becomes diluted to assess the accuracy 

Q = J d^xa^{to)JpQ{to). 


of Eq. (|186|): 


(189) 


























34 


If we assume all dynamical scales are tied to H, we can estimate 


and the diluted charge density at any time t after the initial time to is 

a^{tQ) 


(190) 


JpQ{t) ~ O(iy)(p+(to)<3+(to) 


a^{t) 


(191) 


This translates into an Eq. (186) accuracy of 


accuracy = 0{e 


-3AN\ 


(192) 


where AN > 0 is the efold time from the beginning of the inflation to to the time of horizon exit of a 
given mode. For example, if we want an accuracy of 0(e), we only need the beginning of inflation 
and the time of horizon exit for the longest wavelength mode of phenomenological interest labeled 
by to be separated by the efold number of 


AN ~ —-Ine 
3 


(193) 


before we can set t\ of Eq. (186) to 4 . while assuming for t > 4 . . For e = 0.01, one only 


gives up about 1 efolding. During this one efolding, the (p+(t) decays compared to its value at the 
beginning of inflation (pj^(to). Although its exact trajectory is initial condition dependent, one can 
estimate a lower bound on its decay as long as the initial conditions have < —c+H(p+ which is 
an attractor solution: 

'^scale(^o) 




(194) 


mm 


Since <p+(to) ^ Mp, this sets an upper bound on (p+(4 . ) for the validity of the analytic compu- 


mm 


tation: 


*P+(^fc • ) / 

-min_ < (accuracy)! 

Mp ^ ’ 


1 _ V' 

2 3y 


(195) 


To compute the da correlator using theorem 3, we need the mass of a. We can obtain the mass 
by writing the equation of motion for a, neglecting the small corrections proportional to c_: 


□ la = 0 




2(pi 


-y/gd^(p3 


(P+ 

1 1 
2<P+ y/g 


dp [y/gad^a]+c+H^(p+ = 0 . 


(196) 

(197) 
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The fact that a acquires a time dependent mass is important because that means that the decay of 
the a field due to the mass will stop after a finite time period: i.e. because a is a NG boson, its 
mass will shut off in the vacuum. Hence, we see that in the limit of a/(p+ <C 1, we find 


(□ + c+//^) a 0 

(□ + c+//^) (p+ 0 . 


(198) 

(199) 


There is a remarkable symmetry in this limit because of C/(l)/>g as explained in Eq. (186). The 


mismatch between da and coming from Eq. (196) is suppressed by a/(p+ <C 1. 

Now, once 9 + rolls to the minimum of Eq. ( 184| ), terms proportional to C- dropped in Eq. (197) 
will become active allowing 


(p+ = constant = 

We will call the time when 9 + settles down to this value such that 


( 200 ) 




<P+ 




( 201 ) 


t=tr 


At time tc, a will become massless because it is a NG boson. On the other hand (p+mass does 
not shut off even after 9 + reaches its minimum. Hence, it can be shown (similarly as in ll9^ ) 
that 5(p+ keeps decreasing while da freezes out. Indeed, this decay of the perturbations which 
make it negligible is similar to the reason why one uses a Yukawa interaction to generate ob¬ 
servable isocurvature perturbations in the fermionic isocurvature perturbations of HIOIH : i.e. the 
non-interaction piece has a blue spectrum of n 7. 

Since the axion mass seen in Eq. ( |180 ) is time dependent, corollary 2 is useful to compute the 
isocurvature spectrum. We find 

2 / ^ \ 3-2v{^/^H)-2ek^^+0{ekJgo{m/H) 




( 202 ) 


and the QCD axion fractional density from coherent oscillations can be estimated for Fa -C 10^^ 
GeV as (see e.g. Eq. (14) of HTTH l 

2 / _ /—:-:-\ ’^PT 
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V^^<P+(0 + <P^(^c 
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where 


npT = 1.19 Wa= - 


(204) 


and we have used Q-cdm^^ = 0.12. Here is fixed fidueial wave veetor and is the time when 


that mode leaves the horizon. Furthermore, we ean use Eq. (|171 ) to set 


^ (Plitc) + (pUtc) - 


;(C-+C+) 


y/C-C+ 

eorresponding to the minimum assuming Fa^ H. This gives a CDM fraetion of 


(205) 


npT 




1012 GeV 


(206) 


whieh saturates the relie bound 


Fn ^ 9 


-'l/npj 


(t,J X 1012 GeV. 


'+ 


(207) 


Unlike in the ordinary axion seenario where FI/(InFa) sets the varianee of the effeetive elassieal 
initial eondition for 0+(t^g), here H/{2K^^{tkQ)) <C H/{InFa) sets the varianee.^ Beeause the 


axion is a non-thermal dark matter field after the end of inflation and beeause Eq. (|201[) is expeeted 


to be eontinually satisfied after the end of inflation due to the weak time dependenee of Eq. (200), 
we expeet Eq. (202) to be a good approximation to the final blue isoeurvature speetrum in the 
model of [[H . 

Eor the validity of Eq. ( |202 ) eoming from eorollary 2, the wave veetor k must lie in the range 

<k < kmax (208) 


'^mm 


where and kmax have parameter dependent eonstraints whieh we now diseuss. In praetiee, 
^min ^max should be ehosen to saturate the most stringent of the bounds listed in the eorollary. 
Sinee Buneh-Davies initial eonditions must be set up for the quantum fields after inflation starts, 
we must have > <^{tb)H{th) (where G is the beginning of inflation). Sinee the mode exit time 
h < F (where U is the time that the time dependent mass shifts): 


^min ^ ^rnax 


< X (aeeuraey)- 



2 

3-2v(V?+«) 


(209) 


^ Hence, even with H ~ 10^^ GeV, one can have ^+(fip) ^ Mp and thus tune 0+ as small as 10 ^ without considera¬ 
tions of the variance. Of course, one can even go much smaller than the variance with 10% tuning as well. 
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Figure 1 : Comparison of the result of this paper with the result of fT] for the amplitude of ^ (e g- can 
be taken to be the CMB scale), neglecting small slow-roll parameter corrections of 0{e). The enhancement 


is limited by the condition of Eq. (651. For mjH = 1.497,one can have n 3.8 and an enhancement of 
about 11. 


where we required (p+(f) to be smaller than Mp (“aeeuraey” is defined by Eq. (192) assoeiated 
with the attraetor assumption). If other eonstraints on ^max that we diseuss below allows it, 
^max can realistieally beeome very elose to saturating the upper bound. One ean eonveniently set 
“aeeuraey” equal to the analytie approximation aeeuraey / that one seeks, although in prineiple. 


they ean be set independently. One of the important eonstraints related to / eoming from Eq. (133) 
with k = kmax is 

2 


vmax 






( 210 ) 


Also, beeause one is using a sealing approximation about the fidueial wave veetor k^, there is a 
set of inflationary model dependent validity eonstraints. Given the eomplieated parametrie depen- 
denees of these boundaries, we will summarize below and give an explieit example. 

The main eorreetion to the original isoeurvature result of [[II (i.e. Eq. (21) after aeeounting 
for footnote 5 in [[ 0 ) in the blue speetral index region (i.e. before the break) is the speetral index 
dependent faetor 

A2(k) _ 22KV^")-i|r(v(yc^//)) |2 


Rs{c+) = 




Tt 


( 211 ) 
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where one must keep in mind that v eannot be exaetly zero beeause of Eq. (65).^° This is an 
interesting faetor beeause it ean give a faetor of O (l/(n — 4)^) enhaneement for the blue speetrum 


due to the prefaetor |r(v)f while still satisfying Eq. (65). Eor example, with m/H = 1.497, we 


have n k, 3.8, we obtain an enhaneement faetor of 11. To satisfy Eq. (65), we merely need to 
satisfy 

0.09- > 1 

£k 

for inflationary seenarios with = 10^“^ even if one desires a pereent aeeuraey of / = 0.01. Sueh 
enhaneements are surprising beeause as v —)■ 0, one expeets the quantum fluetuation amplitude 
of a to be smaller beeause its spaeetime-eurvature indueed taehyonie growth is smaller. Indeed, 
one ean easily eheek that {danad^<^nad) indeed is smaller as its effeetive mass eontrolled by c+ is 
inereased. Eor the same parameters, the baekground field a also deereases as c+ is inereased. The 
asymptotie amplitude of {danad^<^nad) assoeiated with the Hankel funetion solution deereases less 
as c+ —)■ 3/2 from below. This is a nontrivial differenee between the quantization indueed Hankel 
funetion versus the homogeneous mode funetion with slow-roll initial eonditions. The eomparison 
of the new result with the old one is given in Eig. 

Next, we will summarize the explieit form of the improved speetrum together with their validity 


eonditions. Using Eq. (205) and Eq. (186), we ean rewrite Eq. (202) as 




R (r ) f 

[iTtCp+itk. 


^ \ 3-2v(^//)-2e*^,+0(eig)go(^/^) 
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( \ 3-2v(v/fT") 
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1/4 

) / 


^min ^k < kmax < k 


( 212 ) 


kr ^ k ^ kp 


''W = W'-9H? 


(213) 


Note that if one wants to compute this in a numerical setting, one must be careful to tune the background field 
boundary condition «(initial) keep the axion background field value at fiducial mode horizon crossing 
fixed for any fixed choice of ko and m. Also, note that as long as Eq. ( |65| l is obeyed, the expansion of the Hankel 
function from which the r(v) arises is a good approximation. 
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(214) 

(215) 

(216) 

(217) 

(218) 


/ '(p+{tko)' 

where A|(fco) ~ 2.4 x 10^^. We have also assumed Fa -C IQi^ GeV for the baekground axion dark 
matter fraetion. Some of the important baekground equation of motion simplifieations allowing 
the analytie treatment eome from 
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>1/4 


kn 

-H<ZFa<Z (p+{tka) < <?>+(?fcHn) ~ ( — 

^Hq 




(219) 
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( 220 ) 

( 221 ) 

( 222 ) 


(c-C+)i/4 

where k/Zp ~ aoHon/l is the wave veetor eorresponding to the observable universe today. The 
veetor kg — a{te)H{te) is the last wave veetor to leave the horizon at the end of inflation and is 
typieally at an unobservably small seale and is very model dependent. This gives the improved 
isoeurvature speetrum (together with Eqs. ( |204| ) and ( 134[ )). The independent variables ean be 
elassified as axion-dependent parameters [Fa, (p+(tz;o), Q+{hQ)F±}^ inflation-dependent parameter 
and approximation seheme dependent parameters {/, ko}. Sinee physies is obviously in¬ 
dependent of different approximation seheme, the physieal parameter spaee of this model is six 
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dimensional. The blue spectral index is controlled by only one parameter c+ while the break in 
the spectrum is determined by kmax 


Eq. (220) gives the approximate condition for classical initial condition tuning of 0+ assum¬ 


ing there being a quantum noise of 0{H/(In)) for the axion field because there is no enhanced 
symmetry that would protect the axion field from tadpole corrections. Although we do not ad¬ 
dress the tadpole issue here by an explicit computation, it is reasonable to expect that the tadpole 
quantum fluctuations can also significantly correct the background equation of motion (again in 
the absence of enhanced symmetries) later as settles to its minimum. In that case, we should 
apply a condition 

H(c^c^y/* (223) 


47tFa^/C-+C+ 
which is stronger than Eq. ( |220| ) because of Eq. (222). 

The isocurvature to adiabatic perturbation ratio is controlled by 


Ajik) 


= PQCD{c+,C-,Fa)9_^_{tkQ)S£kQX 
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'^mm 
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V\/ 2 Ffl ) C++C- \H{tc) F 

where the inflationary scalar spectral index is given by 

ns—\= 2r\y — 6e 


< k.ko < kmax < K 
kr < k < kp 


(224) 


(225) 


which can be used to phenomenologically specify rjy once — 1 and e are fixed. The CDM 
fraction is given by 

tlpT 

^ / 2 Efl,/^ 5 =(c_-t-c+) \ 


(Oa^WaOlitka 


^^1 ^c_c+ ' 


1012 GeV 


< 1 . 


(226) 


Because of Eq. ( |139| ), phenomenologically allowed parameters are in the regime of (0^ 1. If one 

accomplishes this with bringing down Fa, then Eq. (|219[) brings down FI. This in turn brings down 


<P+(tfco) because of Eq. (218). 


Eet us illustrate this formula with concrete parametric choices. We plot A2(k)/A|(ko) in Eig.j^ 
for c+ G {0.2, 1, 1.5} with the rest of the parameters fixed at 


c_ = 0.9, 0+ = 0.04, Fa = 7.9 x lO^'" GeV, (p+itk^) = 8.3 x 10 ^H = 6xl(f GeV, / = 0.2. 

(227) 
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- 1=2 - 3000 

- n=1.1 

. n=1.8 

— - n=2.3 


0.001 0.010 0.100 1 10 100 1000 

kIkQ 

Figure 2: Illustration of the blue isocurvature spectra for several values of the parameter c+ G 
{0.2, 1, 1.5} which map to isocurvature (spectral index, axion dark matter fraction cOa} of {{1.1,4.73 x 
10-4}, {1.8, 4.06 X 10-4}, {2.3, 4.13 X 10 4}}, respectively. The rest of the parameters are fixed as given 
in Eq. ( |227[ ), where in particular, the dark matter fraction is. The gap in the dashed curve in the range 
k/ko G [10,100] occurs as a result of breakdown of the analytic approximation associated with the fact that 
the effective time-dependent mass transition occurs before the modes classicalize. Note that the c+ controls 
both the inflation induced Hubble scale mass for the non-vacuum axion as well as the amplitude of the 
spectra. That is why the three curves do not meet at a point. The dot-dashed CMB curve represents a flat 
spectrum with an amplitude of about 2% of the adiabatic spectrum extending from 1 = 2 — 3000 scale with 
ko = 0.05 Mpc-'. Because of the transfer function suppressing isocurvature power relative to the adiabatic 
power on short length scales (an effect not shown here), the observational constraints on these example 
spectra are weak. 

This set of c+ map to {spectral index, axion dark matter fraction cOa} of {{1.1,4.73 x 
10^"^}, {1.8, 4.06 X IQ-^}, {2.3, 4.13 x IQ-"^}}, respectively. To compare with the approximate 
CMB length scales in the plot, we have fixed ko = 0.05 The gap in the dashed curve 

in the range k/ko E [10,100] occurs as a result of breakdown of the analytic approximation as¬ 
sociated with the fact that the effective time-dependent mass transition occurs before the modes 
classicalize. The actual spectrum in this gap is not addressed by the techniques of this paper. A 

' * Given that this paper is a paper focused on anafytic computation of the spectra, we feave more detaifed numericaf 
data fitting work to the future. 
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- CMB 

— COa=2x^0~^ 
■■■ ai^=6x10-^ 


— - uncorrected a)a=2^x^0 ^ 


Figure 3: Isocurvature spectra for Fa = {6.8 x 10^° GeV, 1.6 x 10^^ GeV} corresponding to the axion dark 
matter fraction of (Oa = {2x 10^^, 6 x 10^^}. The other parameters are fixed at values given in Eq. (2281, 


and the resulting isocurvature spectral index is n = 3.8. The “uncorrected” label refers to the plot that would 


have been given based on the previous literature without Eq. (2111. Note that without the correction, one 
would misidentify an experimental signal of 0 )^ = 2 x 10^^ for that of cOa = 6 x 10^^. The CMB label and 
the gaps in the spectra are explained in Eig. Note that parameters were chosen such as to stay close to 
being observable. It is easy to choose parameters such that this spectrum is unobservable in the foreseeable 
future. 

similar breakdown of the assumptions leads to the termination of the dotted curve. 

Next, we illustrate what happens when the isocurvature spectral index is very steep such that 
the correction factor of Eq. (211) becomes around a factor of 10. In Fig. we plot the isocurva¬ 
ture spectra for Fa = (6.8 x 10^® GeV, 1.6 x 10^^ GeV} corresponding to the axion dark matter 
fraction of C 0 a = {2x 10^^, 6 x 10^^}, with the other parameters fixed at 


c+ = 2.235, c_ = 0.9, e+ = 10“^ (p+{ho) = H = 9 xlO"^ GeV, / = 0.7. 


(228) 


The relatively extreme c+ parametric choice gives an isocurvature spectral index of n = 3.8 and 
the spectral amplitude has an analytic approximation error at the level of around 70%.^^ It is 
remarkable that dark matter fraction Wa as small as those considered in Figs. andcan generate 
potentially observable effects in cosmology. 


12 


A smaller choice of / (corresponding to a smaller approximation error) leads to the spectrum not being computable 
analytically in the interesting region. 
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Another question one might have is how large H ean be in these seenarios, sinee H eon- 
trols the amplitude of the tensor perturbations that may be observable by future experiments 
[|102l4107ll . Unlike in the the ordinary axion situation, the effeetive PQ symmetry breaking 
VEV is mueh larger throughout the eourse of the isoeurvature field evolution. This means that 
///(time dependent PQ symmetry breaking VEV) whieh eontrols the effeetive amplitude of the 
isoeurvature is mueh smaller for observable k region for the same value of Fa whieh eontrols the 
axion dark matter fraetion (Oa- The parametrie tension still arises as we will now see. Combining 
the stabilization eonstraints of and b (i.e. Eqs. (|178 1 and (179)) with the dark matter bound of 


Eq. (226) gives 




/ 2H/h 


npj 


<««<!. 


(229) 


Vl0i2GeV^ 

Another eonstraint eomes from the isoeurvature bound of Eq. ( |224 ), where we restriet to the 
relevant subset of eonstraint 


A 2 (to) 


< Ctk. 


(230) 


Both the dark matter bound of Eq. ( |229| ) and the isoeurvature bound of Eq. ( |230 ) ean be satisfied 
for suffieiently small 0+{ko). To find the minimum 0+(ko) allowed by the other eonstraints, note 


Eq. (221) (eoming from deeoupling axion mixing) and Eq. (223) (eoming from the negleet of 


quantum tadpoles), and Eq. (220) imply 

Henee, there is a minimum c_ for whieh this ean be satisfied: 
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Putting this into Eq. (|231|) gives 
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in(0i) = F (§) = 4 I 2-/3 (y(H/F„)) ‘/3 + 
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(233) 


(234) 


(iy(H/F„))‘'/ 

where the right hand side is valid whenever it is smaller than unity. Now, the isoeurvature bound 
in the form of Eq. (|230|) eombined with Eq. (234) give 


W^2^^PT- 
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The Fa dependent funetion on the left hand side of Eq. (2351 is a monotonically increasing function 
of Fa and it is also a monotonically increasing function of FI. Hence, to maximize FI, we need to 
determine the smallest Fa that is allowed by our constraints. 

One blunt constraint comes from the fact that /i is supposed to be a small angle. Since the small 
angle squared il{Fa/Fl) is a monotonically decreasing function of Fa, one can find the smallest Fa 
allowed by the small angle assumption by solving ii{F\/Fl) = 1: 


F\ 

^ 0.05. 
H 


(236) 


Another constraint comes from combining Eq. (229) with Eq. (234). Defining F 2 = min(Efl) sub¬ 
ject to this constraint, one finds numerically 


Fi 

^« 0.2 


(237) 


which is larger than F\/Fl. One also finds numerically that /z > 0.3 up to the perturbative limit in 
this corner of allowed parametric region. Because the Fa dependence on FI bound is weak, there is 
only a small shift in using F 2 versus F\. In any case, since F 2 is a stronger constraint, we set Fa = F 2 
in Eq. ( |235[ ) and need to solve for the H upper bound. Note that setting Fa = F 2 corresponds to a 
min ( 0 ^) = 0.04 and C- = 0.08c+.^^ 

To solve for the FI upper bound, we still need to set c+. Note although c+ ~ 1.3795 gives the 


smallest Eq. (211) that shows up in Eq. (235), we must still check the constraints coming from 
Eq. (219|). Eq. (178) part of the constraint in this parametric corner can be written explicitly as 


(c+c- 


1I/4 


h 


-H^Fa 


(2At(0.2))i/4vT^ 


h 


< 0.2 


(238) 


where one sees the explicit H independence because of Eq. ( |237| ). Another can be written as 

<P+(^/to) 


Ma 






Einally, perturbativity requires 


J^Hq 
h < a/4^. 


(239) 


(240) 


Note that this does not mean that the smallest 0^ for all parts of the parameter space is 0.04. It is only when Fa is 
minimized subject to the constraints discussed, do we have this minimum on 0^ . 
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. f=0.1 a=3x10-2 

- f=0.1 a=3x10-^ 

- f=0.9 a=3x10-2 


Figure 4: Maximum H as a function of the isocurvature spectral index n. On the right axis, the tensor-to- 
scalar power spectra ratio r is plotted. The break feature for small n in the / = 0.9 case corresponds to the 
situation where Eq. (|239l becomes more important than the constraints coming from Eqs. (2421 and (2431 


These three eonstraints constrain {h, c+, <p+(tfeo)} for a fixed non-physics parameters such as ko 
and Eqs. (2^) and (2^) put a bound of 


c+(for max//)< 1.78 (241) 

corresponding to n < 2.6. In this c+ range, the c+ dependence of Eq. ( |235| ) is weak. Erom 


Eq. (2351, we thus conclude that the maximum H that is allowed by the present scenario is 


(242) 


Next, we still need to impose the constraints of Eq. ( |218 1 from the quantization. This constraint 
requires H to be large and (p+/Mp to be small. We can thus denote this constraint as a lower bound 
on H : 

0-82 ^Q_3^-|-rv+o(e)go(m/^) / 9+{%)/Mp 


f 


X 


8.7 X 10 “ GeV<H. 


(243) 


The choice of ko does contain phenomenological information in where the data constraints lie because one usually 
wants to choose ko where the data is accurate. 
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For any fixed c+, both Eqs. (242) and (2431 can be satisfied if (p+itko)/Mp is small enough. With 


the right hand side of Eq. ( 242| ) set equal to the left hand side of Eq. ( |243| ), we can find another 
constraint on maximum <p+(t/to) /^p similar to that coming from Eq. (239). This latter constraint 


usually is more important than Eq. (239) for a blue spectral index. We thus arrive at the maximum 
H shown in Eig.|^for the blue spectral index scenario considered in this section. The correspond¬ 
ing tensor-to-scalar power spectra ratio r is also shown on the right axis of the same figure. It 
is clear that inflationary scenarios consistent with a very blue isocurvature spectra (e.g. n > 2.0) 
will not generate tensor spectra that is observable in the near future. Conversely, a tensor-to-scalar 
ratio at the level of r = 0(10^^) will disfavor this class of models. 

Eet us now summarize the how one obtained the upper bound on H. One minimizes the angle 
0 +(^/to) consistent with the validity of the classical equations of the motion and subject to the 
decoupling constraints. Putting this into the isocurvature formula and minimizing the isocurvature 


to curvature ratio by varying Fa/H subject to dark matter abundance constraints lead to Eq. (242) 
bound on FI. The decoupling constraints are not particularly fundamental, but the study of that 
region is beyond the scope of this paper. Hence, even higher values of the tensor-to-scalar ratio 
may be valid with blue isocurvature perturbations, but the phenomenological signatures will be 
more complicated than the simple situation presented here. 


3.2. Do Dressing Effects Give a Lower Bound on the Blue Isocurvature Spectrum? 


Erom the definition of Sy in Eq. (17), we would naively expect 




-I- cross terms 


(244) 


where the second term arises from the “dressing” effect of the isocurvature coming from the fact 
that spectator isocurvature is a contrast between the isocurvature field and the adiabatic field. Eor 
scale invariant isocurvature spectra (such as massless axions), we have Zo = 0 which means that 
this second term coefficient is negligible. However, for a blue spectrum, this coefficient is of order 
unity. Hence, one would naively expect 


{S^S^) > (CO naive operator product expectation for blue spectra 


(245) 
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independently of /Xo amplitude.Henee, a little puzzle arises whether one ean 

conelude 

A? 

(246) 


^ > 1 
A 2 ~ 


for a blue isoeurvature speetrum in the limit H 0 with miH = 0(1) and Xoi^ko) fixed. 
Theorem 3 trivially states that this lower bound with Xoi%) fixed does not exist sinee 


-4- ^ 1.3 X 10^ 


x22''x |r(v)|2(^ 


10 


Mr, 


(247) 


whieh vanishes as the inflaton slow-roll parameter e —)■ 0 with fixed. That means that 


(CO 


^0 


(248) 


for blue speetra is possible if Xoiho) can be fixed eontrary to the naive expeetation from 
Eq. (245).^^ On the other hand, consistent quantization does give a lower bound shown in 


Eq. (139). However, Eq. (139) arises because Xoiho)/Mp cannot be fixed because of the quan¬ 


tization approximation used: i.e. the reason for Eq. (139) is different from the naive operator 
multiplication analysis. Indeed, (ko) shown in Eq. (63) is independent of (up to the inequal¬ 


ities implied by Eq. (69)) while A^(fco) shown in Eq. (14) manifestly does depend on 

If one looks at the details of the proof to try to understand where the naive operator multiplica¬ 
tion analysis goes wrong, one notes that there is a secular growth effect in which the isoeurvature 
field Sx^^^ develops an adiabatic piece due to a secular superhorizon source effect. Eor exam¬ 
ple, in the spatially flat gauge, the gravitational interaction transmitted adiabatic piece is given by 


Eq. (93). The fact that this comes from gravitational physics can be understood from the fact that 
the off-diagonal mass squared in Eq. ( |8^ are Planck suppressed, and these terms are responsible 
for Eq. ( [93] ). This can be interpreted as the effect of gravity imprinting dominant energy inhomo¬ 
geneity information onto the subdominant isoeurvature field. Hence, even if 5x^^'^ at the quantum 
fluctuation level does not have ^ correlation information^^, it will on far superhorizon scales look 


This does not by itself give a bound on the total isoeurvature perturbations which depend on the dark matter fraction 

This limit e —^ 0 can be fraught with strong coupling issues in the density perturbation computation formalism. We 
can neglect these issues and can take this limit formally since the point is that it decreases towards zero and not 
about the absolute magnitude. 

In the proof of theorem 3, the spatially flat subhorizon modes are essentially decoupled from the inflaton modes 
by Xo/Mp. This allows one to determine the Bunch-Davies state quantum correlator independently of the sourced 
mixing with the inflaton in the subhorizon region. 
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like a mixture of adiabatie and nonadiabatie field, in preeisely the eombination to eliminate the ^ 
dependent pieees in (S^S^). 


4. SUMMARY AND CONCLUSION 


In this paper, we have presented three theorems and related eorollaries eoneeming blue spee- 
tra produeed by linear speetator isoeurvature fields that give rise to CDM-photon isoeurvature 
perturbations. Theorem 1 defines a superhorizon eonserved quantity for systems possessing an 
approximate symmetry of ^ V^{5x)Sx- The merit of this theorem eompared to previ¬ 

ous diseussions of this topie in the literature is its ability to go beyond the end of inflation and 
the reheating proeess. Theorem 2 deseribes under what averaging eonditions that fluid quantities 
behave as dXnadIXo- This seeond theorem merely restates what is known in the literature (see 
e.g. [I^ |92l) in the eontext of eurrent theorems. Theorem 3 deseribes the eomputation of the 
quantum isoeurvature perturbations. The merit of theorem 3 eompared to the previous diseussion 
in the literature is the explieit eanonieal quantization in the presenee of linearized gravitational 


eonstraints. The validity regime of this theorem imposes a nontrivial eonstraint of Eq. (138). If 
this eondition is violated, the amplitude of quantization is expeeted to be more eomplieated than 
the simple analytie treatment presented here. 

In See. 3.1[ we have applied the theorems to the work of [HI and improved their eomputation. 
In the proeess, we have uneovered a eonserved Noether eurrent assoeiated with U{l)pQ that is 
leading to the traeking of the axion field with the radial direetion field. The final speetral formula 


is given in Eq. (212). The general magnitude eomparison of the isoeurvature blue part of the 
speetrum is shown in Eig. The speetral break features and the validity of analytie eomputations 
were explored in Eigs.|^and[^ The maximum tensor-to-sealar ratio for whieh this simple seenario 
remains valid is shown in Eig. 

In See. 3.2[ we have applied the theorems to explain how naive operator produet estimates for 
the isoeurvature eorrelator amplitude lower bound fails. The main physies is that the speetator field 
attains the inhomogeneities assoeiated with the inflaton through its gravitational eoupling. Erom 
a perturbation theory perspeetive, this inhomogeneity is attained through a seeular effeet whieh 
would naively be dropped from the eonsideration of perturbative expansion eoeffieient alone. Erom 
a physieal perspeetive, the speetator field whieh undergoes no appreeiable quantum fluetuations 
by themselves still attains an inhomogeneity that looks like the inflaton’s inhomogeneities. Inter- 
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estingly, we do uncover in this paper an isocurvature correlator amplitude lower bound Eq. ( |139[ ) 
whose phenomenological validity requires the dark matter fraction co^ to be much smaller than 
unity. If this lower bound is violated, the quantization of the isocurvature perturbations do not take 
on the simple form presented in this paper. 

There are many possible future extensions of this work. Regarding the general applicability of 
the theorem, it would be interesting to find interaction strength boundaries for classes of models 
for which the linear spectator behavior of the isocurvature perturbation survives. Regarding the 
scenario of m, we have explicitly laid out where the analytic computation fails near the break 
region of the spectra. Although we would naively expect that either side of the break region to be 
smoothly connected, we would also naively expect features to exist in that region. Numerical in¬ 
vestigations of the features may be interesting for discovery potential. Other obvious future inves¬ 
tigation possibilities include improving our understanding of the experimental discovery prospects 
of the blue spectral isocurvature perturbations. 
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Appendix A: Particular Solution In the Subhorizon Region 


Consider the approximate equation of motion object for coming from Eq. (781 



where we have approximated 



(A2) 


Without taking the superhorizon limit, consider the following particular solution ansatz: 




(Po 


(A3) 
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The equation of motion objeet E becomes 




Xo (Po 


S^ + WSq) j 


Zo Mj Hxo (po 


V^^o) yXoVipicPo) ^ ^KicPo)? 


Xo(po 


‘Po 


Next, the usual slow-roll approximation gives 

(po 


-6H 


which implies 


Zo 


IT ^ ^+MI 2 -+ (3t]v - 6 e)//^ 


Note that Eq. (A7) assumes a small slow-roll factor analogous to £ but for the Zo field 


call these expansion parameters 


m2 , _ 1 Zo 


'‘lx = % = — TTT TTV < 1 


547/4^2 


whose motivation is detailed in Sec.|^ Explicitly, one can show 

+ = 0(£^^^)+0(^/£r|v) + 0(^/£r|^) 

Xo <Po L 

Next. we know 


H 


which yields 




W ^ ^ + {3riv - eey 


Einally, we also write down the equation of motion for as 

(f5y 


^^2 +3//<9,5(p(^^)+ (^^ + (3T]v-6£)//2j 5(p(.^/)+;„2^5^(./) 


and note that 


2 ^0X0 ^{sf) _ ^2 ^ V^ZO ^ (sf) _ ^ ( Xo 


nr ' = m^sgnffln 

MjH ^ ^ ^ M, 




:.ri) 


We thus conclude 


E = 0 + O(^^]+O (£”>i) + O(Jiy) 


(A4) 

yi ( Po ). 

(A5) 

(A6) 

(A7) 

(A8) 
. We will 
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(A12) 

(A13) 

(A14) 
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which means that Eq. (A3) solves the equation of motion even in the subhorizon region in the 
leading approximation. 

Finally, note that the partieular solution itself is of order 




<Po 


MpJ 


(A16) 


which means that sinee we are dropping the same order in Eq. (|A15 1, one might naively think 


there is no eontent in the solution. However, note that on the left hand side of Eq. (All, there is a 
factor of {k/a)^ which makes this term non-negligible. The point of this seetion was to show that 
sueh unsuppressed terms are all eaneeled by the slow-roll equations of motion to leading order in 
slow-roll expansion. 


Appendix B: Mixture 


Result of theorem 3 is applieable only when the dark matter x can be made to be totality of 
dark matter. Suppose 

^Pcdrn Spx + Spy 

-=-^ 

Pcdm PX + PY 

where X ean be the x partiele and Y is the rest of the eold dark matter. We ean rewrite this as 

Spcdm _ Px^x + pySy 
Pcdm Px + pY 

= (Oxdx + CJydy. 


(B2) 

(B3) 


where 


sueh that C0x + 0>y 


1. Henee 


(Ox = 


Px 

PX+PY 


(Oy = 


PY 

Px+pY 


_ ^Pcdm 3 Spy 

Pcdm 4 py 

3 

= (Ox Sx + (Oy Sy — -Sy 


(B4) 


(B5) 

(B6) 


If 


then 



Ss = cOx 



(OxSsx 


(B7) 

(B8) 


where Ssx is the isoeurvature in X component. Hence, with mixing, the isoeurvature is diluted by 
a faetor (Ox- 
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Appendix C: Slow-Roll of x 


In this section, we motivate a slow-roll expansion parameter for x field. Consider 

Zo + 3//;i:o + Vj(Zo) =0 

where H is an externally determined time dependent function. Divide through by HM^. 

Zo ^ gZo ^ _Q 


HMl Mj HMj 


Define a formal hierarchy 


Zo 

HMj 




Zo ^;^(Zo) 


m2’ HMj 


and expand Xo/^p about zero using a formal perturbation parameter A: 


f^ = n(X) 


= 0{X)} 


Using this expansion, construct a trial solution 


Zo = 


3H{t) 


X + vi A^. 


(Cl) 


(C2) 


(C3) 


(C4) 


(C5) 


Put the trial solution Eq. ( |C5[ ) into Eq. ( |C1[ ). Eor this endeavor, we need to evaluate Zo to second 
order in A: 


Zo = 




3H{t) 


3//2 


V^'(Zo(t)) V^(Zo(0) ^2 ^ + D(A3) 


(C6) 

(C7) 


3H{t) 3H{t) ' 3//2 

To finish expanding the right hand side, we need to evaluate H. According to the usual slow-roll, 
we have 


H = -£H\ 


(C8) 


Hence, we conclude 


Zo = 
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9H2 X 


V'MK(Xo)>--- + ^^>1 + O(A’) 


(C9) 


Put the trial solution Eq. (C5) into Eq. (C2), account for the second term of Eq. (C4), treat 
eX = D(A^) in the formal counting, and collect 0{X) and D(A^): 
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This can be solved 

vi = f dt 


KUo) 


i3 


54//4M2 




Hence, we conclude the fractional correction to ^ is 


vi 


K(xo{i)) 

3H{t) 
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1 
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vlcixoii)) 3H 5AH^Ml 

\ m) ^ 


This motivates us to define the following slow-roll parameters: 

^ X 3^ X (^o) = [xq] = vx 


3H 


Vxixo) 3H 9//2 ^ 

1 1 

X -— X 
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(C14) 


Obviously, this is not unique, and other slow-roll definitions exist. See e.g. 
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